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v ABSTRACT  L — “ 

The  coherence  between  two  points  of  the  received  acoustic  wave- 
form after  reflection  from  a time  varying  random  surface  is  evaluated  for  the 
far  field  case.  A Neumann-Pierson  spectrum  and  an  isotropic  sea  is 
considered.  For  the  low  roughness  case,  the  coherence  is  computed  for 
wind  speeds  from  2 to  10  knots.  For  the  specular  and  non- specular  direction, 
the  coherent  region  is  on  the  order  of  the  pattern  resulting  from  the 
insonlfied  surface  area. 
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INTRODUCTION 


A very  important  problem  in  underwater  sound  detection  is  the  com- 
putation of  array  gain.  The  earliest  calculation  of  array  gain  assumed 
that  the  signal  was  perfectly  coherent  and  that  the  noise  between  re- 
ceivers was  perfectly  incoherent.  Thus,  the  signal  added  as  N^  and  the 
noise  as  N,  where  N is  the  number  of  elements  in  the  array.  The  array 
gain  is  10  log  N for  this  idealized  situation.  Later  work  assumed  that 
the  noise  was  isotropic,  that  is  the  same  in  all  directions.  For  this 
case,  the  noise  was  incoherent  at  half  wave  length  spacing  and  almost 
incoherent  when  spacing  between  elements  is  greater  than  2 wavelengths. 
This  model  was  later  extended  to  the  case  of  directional  noise  which 
showed  that  the  array  gain  was  very  dependent  on  steering  direction. 


In  comparison  with  the  noise  models,  relatively  little  work  has 
been  done  on  signal  coherence.  A notable  exception  to  this  is  the  ex- 
cellent work  by  Parkins1.  In  this  memorandum  a generalization  and 
correction  of  the  work  by  Parkins  on  small  roughness  is  made.  In  ad- 
dition, plots  of  coherence  are  provided.  The  purpose  of  this  memoran- 
dum is  to  provide  a model  and  computations  for  some  typical  cases,  for 
the  coherence  of  a signal  reflected  from  the  ocean  surface. 


REFLECTED  PRESSURE 


We  will  consider  a single  frequency  waveform,  exp  (12irft)  insoni 
fying  a finite  area  of  ocean  surface.  The  reflected  pressure  in  the 
far  field  is 


An  equation  similar  to  equation  #1  was  utilized  by  Eckart  . A dis- 
cussion of  equation  #1  and  the  corresponding  assumptions  Involved  for 
the  rough  surface  case  are  given  by  Nuttall  and  Cron 3.  in  equation 
#1,  Pi(t)  is  the  reflected  complex  pressure  at  the  ith  point  at  time 
t.  B is  a geometric  factor,  X is  the  acoustic  wavelength,  R^  is  the 
distance  from  surface  origin  to  the  1th  point,  k is  the  wave  number. 
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p (x.,y^)  Is  the  insonlfled  pressure  on  the  surface  at  y^;  a^, 
bA  and  c^  are  the  sum  of  the  direction  cosines  of  the  incident  and 
reflected  pressures  for  the  x,  y and  z directions,  respectively. 

T(*i.  y±>  is  the  surface  height  at  the  point  x*,  y^  at  time  t. 
The  incident  mathematical  signal,  exp  (12wft),  has  positive /frequen- 
cies only.  This  signal  is  modulated  by  the  time  varying  ocean  sur- 
face. Since  this  time  variation  is  slow  in  comparison  to  the  varia- 
tion of  exp  (i2*ft),  then  Pi<t)  is  narrow  band.  This  P^t)  contains 
only  positive  frequencies.  P^(t)  is  an  analytic  signal.  The  real 
part  of  P,(t)  is  the  actual  pressure  and  the  imaginary  part  of  PA<t) 
is  the  Hilbert  transform  of  the  real  part. 


MUTUAL  AND  COMPLEX  DEGREE  OF  COHERENCE 


He  now  consider  two  points  of  the  pressure  field.  (See  Figure  1) 


as  defined  by  Born  and  Wolf  is 


The  mutual  coherence  function 


7Z.M  - 


Where  * is  the  complex  conjugate  operator 
< > is  the  ensemble  average 
'T'  is  the  time  delay  between  points  1 and  2 

The  complex  degree  of  coherence  is  defined  as 


By  the  use  of  Schwarts's  Inequality*  it  can  be  shown  that 


o 4 Y«tnl£  i 

The  upper  limit  corresponds  to  a perfect  coherence  of  the  pressures  at 
points  1 and  2,  whereas  the  lower  limit  signifies  perfect  incoherence 
between  the  points  at  1 and  2,  In  signal  processing  terminology, 

is  called  the  complex  correlation  of  Afe) and  +>*.[£}  • It  can 
be  shown  that  • 

TZ.  M - 2 jT  J?aO)  + L 7?/  (r-)J  * V 

Md  (>)  , <^e 

Where  Re  and  Im  signify  the  real  and  imaginary  parts,  respectively. 
Thus,  from  equation  4,  the  correlation  of  two  signals  at  points  1 
and  2 may  be  obtained,  since  . \ 

■ . fyi  = -t  ft*  (r-M) 

MUTUAL  COHERENCE  FUNCTION  FOR  A TIME  VARYING  SURFACE 
From  equation  #1  and  referring  to  Fig.  1, 

Afcl  - expfai q/a  1 n,  <e, 

^ «o  v 

+ C,};  ~ c »r,)J 

In  order  to  keep  the  symbols  to  a minimum,  we  will  set  ^-O  i and 

We  now  take  the  ensemble  average  of  SL  (£)  A and  we  will  drop 

the  t's.  Then  7 


4 
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kfff/jx,  jj,  jxK  jgl  fic*„  u 4lx“ *6 

. expfi-i fa, -a,*.-*- i, {/, - fc^)J ^expCU  (c, x - c,x)J > 


He  will  assume  thac  the  surface  height  is  Gaussian  distributed.  Kins-* 
man*’  cites  experimental  data  showing  that  the  surface  height  is  very 
close  to  a Gaussian  distribution. 

oO 

^ e x. p £ CctX -c*Ti.)] />  e*rD4Cc,rt-cxTj)]fifr,,T,) 

. ~°o 

and  for  a joint  Gaussian  Jb  (Tt/  T^)  this  becomes 

<exP cin.)J>  = exr£-  £?fYc' ■'-cS-zr.  CieSj  *7 

where 

e - c-c  t,-  <=*-) 

substituting  7 into  equation  6,  we  obtain 

expZiJkCa-.x,- ^K  + t,  y, - yj] expCr  w, Cte] 


SMALL  ROUGHNESS  CASE 

In  this  section,  we  consider  the  small  acoustic  roughness  case. 

For  this  case  we  can  expand  the  term  exp(”4sV,vC(  ctP_J  in  equation  #8 
as 

€xp[^v*-c,  c_<?7  ^/+ivLc,ct  e 

Since  ^ ~ , this  expansion  is  sufficiently  accurate  ^ *\. 

^ I . Substituting  equation  #9  into  equation  #8,  the  first  term  of 
equation  #9  results  in  a coherent  component  and  the  second  term  re- 
sults in  an  incoherent  component. 

The  first  term  called  the  coherent  component  is 


5 
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, K ^lr'Cc?+c?)J 

* fffj Jx,  Jy,  jf*  (*•>  I/O  $(•**>  y^) 

* jo  Jk ~ y>~  y^l^J 


&/o 


The  Incoherent  term  is 


<A  A%  - < e*  i‘£-  ^c'  + * " 

•J5yiv‘c,c.f  /fo,  *.)#*«&>  exp['iAct>.,x,-«.t)cl+u,y,-iHi/4J 


The  coherent  component  in  equation  #10  can  be  integrated  £or  some  values 
of  In  this  study,  it  will  be  integrated  for  Gaussian  insoni- 

fication.  The  incoherent  component  is  more  difficult. 

Let  us  now  consider  equation  #11.  Ue  will  assume  that  the  proper- 
ties of  the  surface  depend  only  on  the  difference  of  the  coordinates. 

Let  <J  - y.~  z, 

i ' c 

V = %-ih. 

the.  _/^  = K e*H?  -Q?'’(cl'-+cZ'2  -&  c>-  CC°j  K °) 

oo 

• ffff  do dV + 

Let  ^ 

then 

<A  At^  - X cxpf  ^^,'fcrtJiVV,  ct  Pfy,  ^ °) 
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*10. 


• fff/d4'l$Ao-lv  /<y+J£,P+-£)jfu-£  p-K.\ 

— «>o  1 


We  now  assume  that  the  effective  extents  on  the  surface  of  the  Incident 
Illumination  Jp  are  much  larger  than  the  distances  at  which  the  sur- 
face heights  are  statistically  dependent  on  each  other.  That  is*  the  ' 
correlation  distance  is  much  less  than  the  insonification  distance. 
Using  this  fact,  we  may  approximate  equation  #12  by 

- K **p[-  4?^  *'s 

ecy,  o&xpfuC(^)  " +(^)  vjj 

•fSMi?  |3)  ex  1-- vw  + C t.  - *0  $l] 


Equation  #13  is  a generalization  and  correction  of  Parkin's  work. 
Equation  #13  is  symmetric  with  respect  to  points  1 and  2,  whereas 
Parkin? s results  are  non-symmetric.  For  the  case  of  Gaussian  insonifi- 
cation, Parkin's  equations  result  in  Incorrect  coherence  values. 

Thus  for  the  Incoherent  component,  a four  fold  integral  has  been 
approximated  by  the  product  of  two  double  integrals.  Equations  #10 
and  #13  are  the  general  equations  for  the  coherent  and  incoherent 
components. 


GAUSSIAN  INSONIFICATION  AND  NEUMANN-PIERSON  SPECTRUM 

We  will  now  obtain  the  equations  for  the  special  cases  of  Gaussian 
surface  insonification  and  a Neumann-Plerson  height  spectrum.  Let  the 
surface  Gaussian  insonification  be 


Thus  the  insonified  pressure  is  I at  the  origin  and  falls  to  a value 
of  *c  « distance  of  L units  from  the  origin. 
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Let  us  consider  the  integral 

<=o 


T/  ^ a,*,) 


Completing  squares  or  using  the  analogy  of  the  characteristic  function 
of  a Gaussian  distribution,  we  obtain 


-X,  - fr^Z?  ex  p 


Thus  for  equation  #10,  we  have 


/4, 4f>  „ = kV=>  rr  L^exp[-  41^  (<+?  + 0.1  + k.N-  £)]  *,£ 

'COM  a J 


Let  us  now  consider  the  second  double  integral  in  equation  #13.  This 
integral  is 


oo 


~ f/M  JL  z(ott  (?)  €xp4'cJ^J^.t~0.z)o(  + (b‘-  i>2.)  41 

exp[- 


Again 


Using  the  results  as  given  in  equation  #14,  we  obtain 


Tz^rrc-  ex  e*p£L  Axg,  - ^ =*/£ 


Let  us  now  consider  the  first  double  integral  of  equation  #13. 


r3  =ffe(o^q-)  e»p[q  " + (*■  *•  -i"  Av 


For  brevity  purposes,  let 


X"  - c^‘  +“  ^a.  In 

q,  - j 


r - ^ 


We  now  express  surface  spatial  correlation  0Cujv/O)  terms  of  the 
directional  wave  spectrum  AxCkx/  k.y.')  . and  Ky  are  the  surface 
wave  numbers  in  rectangular  coordinates.  From  Kinsman^ 


I 

1 I 
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Thus  oo 


- ffdu  Jv  expj^c  A ("a.  o 4-  ^ y)J 

•ff  4-h^  h Ax(kz.  k#')  c os  C kx.°  + v) 

- OO  J 0 


u -m<  v')--  + *’e«p&iCk'‘c,+tyv)J 


cosCk*  u + Kyy)-  ^ 

then 

JT3-f/J  **  1 kw  A'C*x,Ky) 

or  * — 

•/  Lexp{c[cJkH-^)^  *?)$ 

— -v> 

4-  exp^fii  -kx)^  + (A  l- kj)  V^J^J 

Integrating  over  u and  v and  using  the  relation 

CO 

/exp6'^)^  - a ?r  «f60 

v cO 

where  Is  the  Kronecker  delta,  we  obtain 

ar3-=  ffJkxJkj  AXk*,k,)(zTr)'lCf(s-4  + *x')S(^+k?'> 
+ SY-lT-Ky'jJ 

Integrating  on  1^  and  Ky,  we  obtain 

t3  - ■*" 

In  order  to  obtain  some  numerical  values  of  coherence,  we  will  assume 
an  isotropic  sea,  that  is  waves  equally  likely  in  all  directions.  To 


r 
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obtain  this*  let  us  first  state  the  relation  between  the  directional 
wave  spectrum  in  rectangular  and  polar  coordinates. 


Then , (See  Kinsman). 

(Kx  + Kf)3/* 

Then  

~r  _ fF  


For  an  isotropic  surface  ' A*  ^ 


J,: 


■rr  f AUfjxc^^^ 

mrl  L ~W TT^P^  _ 


it/9 


3r^rc  (&  * £ ) * j 

Thus,  we  have  evaluated  I3  for  a general  directional  wave  spectrum 
(equation  #17)  and  for  an  isotropic  sea  (equation  #18). 


For  the  Neumann-Pierson  spectrum 

where  g is  the  acceleration  of  gravity 
s is  the  wind  speed.' 
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EQUATIONS  FOR  COMPUTATION 

For  Che  reader's  convenience,  let  us  now  collect  Che  equations 
needed  for  computation.  As  stated  previously,  for  Gaussian  insonifica- 
tion  and  an  isotropic  sea  with  a Neum&nn-Pierson  spectrum  at  a point 

< A A*4W-  k'  n'tX  exi”J}  ‘fj exrjZ 

. i r 4 v*e, c,  , 

n-J-V*  L -1 


r*c,c* 


where 


K,  - — — <sxrC£*r£>')exrp4{%,-#t.jl  ejcpP  4L£*Xc,'-±ci) 
X-R>  ^ J L *-  J 


The  other  terms  have  been  defined  previously 

t f*  exp  6^.)  . 

Let  'TXo  and  consider  points  1 and  2 to  lie  on  a circle  with  the 
center  at  the  origin,  so  that  Rj-Rj. 

In  the  cgs  system  g*>980. 665  cm.  If  s is  expressed  in  knots,  then 
to  change  to  cm/sec.,  we  must  multiply  by  185,200/3600.  In  equation 
#22,  c- 30, 500.  If  the  wind  speed  is  specified,  then  the  mean  square 
height  of  the  surface  may  be  obtained.  This  is 


suffer 


Thus  for  the  Neumann-Pierson  spectrum,  the  mean  square  height  is 
proportional  to  the  5th  power  of  the  wind  speed.  . It  stamld  be  noted 
that  although  we  have  separated  the  equation  for  v*.  into  coherent 
and  incoherent  components,  the  computer  program  will  add  the  two  fac- 
tors together  with  the  corresponding  constants.  To  find  the  complex 


11 
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degree  of  coherence.  It  is  necessary  to  evaluate^/*1  -A'  ^ and 

-A#*^  is  obtained  by  setting  S2“al»  b2*b^  ana  C2*c^.  For  example 
from  equation  #20,  we  obtain 

''  kl* rrLX)  (a < a >>^Sj 

The  computer  program  for  the  computation  of  the  complex  degree  of  co- 
herence is  given  in  Appendix  A. 


RESULTS 

For  our  computations,  we  have  chosen  a frequency  of  400  Hr  and  an 
angle  of  incidence  of  45°.  The  incident  and  reflected  ray  ar*  coplanar 
with  the  normal  to  the  surface  at  the  origin.  One  of  the  reflected 
directions  was  fixed  at  9*45°.  (See  Figure  #2). 


Fig.  2 

The  other  direction  was  varied  in  .1°  increments  from  41°  to  49°.  The 
wind  speed  was  held  fixed  for  each  computer  run.  In  figure  #3,  3 
computer  runs  are  shown  for  wind  speeds  of  5,  7.5  and  10  knots.  As 
stated  previously,  9 represents  the  direction  of  the  second  reflected 
ray.  Since  the  first  ray  is  fixed  at  45°,  when  the  second  ray  is  at 
45  , a complex  degree  of  coherence  of  1 is  obtained,  by  definition. 

All  3 curves  are  thus  correct  at  9-45°.  All  3 curves  show  perfect  co- 
herence for  9 close  to  45°  and  all  3 curves  are  zero  when  9 is  *2°  from 
45°.  Physically,  the  ray  in  the  specular  direction  has  a coherent 
component  which  is  high  at  low  wind  speeds.  The  coherent  component  de- 
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creases  as  9 goes  further  away  from  the  specular  direction.  The  co- 
herence between  the  two  rays  is  due  to  the  coherence  between  the  two 
coherent  components.  However,  the  coherent  component  is  almost  0 
when  the  direction  is  outside  of  the  beamwldth  (BW)  of  the  insonified 
region.  (See  Nuttall  and  Cron^,  Equation  #7)  For  an  1>32  A » such 
as  chosen  in  these  computations, 

BW  ^ A degrees  - A_  - /.  7a° 

l.  rr  ttk  * '*  '7 


BW  ~~  Su 
u 


degrees 


* /.  7<f 


o 

Thus  the  BW  is  on  the  order  of  2 . 


As  L is  increased,  the  coherent  region  decreases.  It  should  be 
noted  that  in  the  far field,  a coherence  on  the  order  of  a few  degrees 
may  represent  a large  linear  region  of  coherence  along  an  array.  S 
knots  represents  a smoother  surface  than  7.5  or  10  knots  and  thus  has 
a larger  region  of  coherency. 


Some  of  the  parameters  associated  with  the  case  shown  in  Fig*  3 
are  tabulated  in  Table  1. 


Speed  in 
Knots 

<r 

in  cms 

§4^ 

2 

.27 

.0062 

.000039 

5 

2.6 

.0617 

.0038 

7.5 

7.28 

.1699 

.0289 

10 

14.95 

.3488 

.1216 

15 

41.21 

.9611 

.9503 

Note  that  the  wind  speed  of 
and  therefore  does  not  fall 

15  knots  has  a that  is  close  to  / 

into  the  category  of  small  roughness. 

For  a wind  speed  of  2 knots,  the  surface  is  almost  like  aosmooth 
mirror  and  computation  shows  a coherence  of  / from  41  to  49°. 


It  is  probably  perfectly  coherent  in  all  regions. 


Fig.  4 represents  the  case  of  one  position  fixed  at  the  non- 
specular  direction  of  55°.  For  this  case  it  was  found  that  as  the 
wind  speed  changed  from  2 to  10  knots,  the  curves  did  not  change  much. 


\ 
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The  coherency  drops  much  fester  for  this  case  then  for  the  specular 
direction.  Fig.  5,  represents  the  case  of  one  position  being  held  at 
65°.  The  difference  between  the  2 knot  case  and  the  10  knot  case  is 
hard  to  distinguish  on  the  plot.  An  investigation  of  the  parameters 
shows  that  both  curves  are  dependent  on  the  tails  of  the  surface  wave 
height  spectrum.  However,  the  tails  of  the  wave  height  spectrum  are 
very  close  to  one  another  at  high  surface  frequencies. 


CONCLUSION 


A model  and  equations  have  been  obtained  for  the  spatial  coherence 
in  the  far-field  for  a signal  reflected  from  a time  varying  random  sur- 
face. Computations  and  plots  of  values  for  specific  cases  have  been 
presented.  Only  the  low  roughness  case  has  been  presented. 


Plans  are  being  made  to  obtain  the  incoherent  component  of  spa- 
tial coherence  for  arbitrary  roughness.  The  Fresnel  region  case  is 
also  being  considered.  In  addition  to  this,  experimental  data  of 
elements  of  an  array  have  been  obtained  and  will  be  analyzed  in  the 
near  future. 
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COMPUTER  PROGRAM 


31 

55 


PARAMETER  NP=O0O » NPl=NP+ a , NTsNP+1 » NUNP+2 • N2=up*3 
REAL  Xl(N2)»Yi(N2),8UFFt*<10000) 

implicit  ooucjle  precision (a-h, o-z>  

DIMENSION  CO(NPI) 

format  (5o6*i*o9«4rD5«o)  _ 

FORMAT ( IX* 6016,8/6016. 8/ 0DI6. 8//) 

DEFINE  F4(X*Y)=SQRT (SORT {X*X+Y*Y) ) 

SpMT=50t9D0 

Sp=bi,D0  T..  

CALL  PL0TS(BUFFER»10000»o) 

call  plot(0»o*ioo)  __  

call  plot <o.»o.*-3) 

CCM=30.4bD0  


CFM=3o500.DO 

0=980.66500  

C=500o.DO 

C=C*CCM  

F=400.0o 

Pi=3. 141592653589793*400  


AM0A=C/F 


WN=2.*PI/AmDA 

WN2=MU*WN 


CF=P1/180  *00 

30  READ ( 3 * 31 * ENO=32 ) T 1 * Hi * Tsl » HS1 * HS2 * SK * PL 

pl=pl*c/p 

PL2=PL*PL 

..  6NF=2.*PI*PL2  

GNF=GNF*ttNF 

_ S=SK*185200.00/3600#.00_:  . 

S2=S*S 

SC=S/l2,*G)  

SC2=SC*SC 

SC4=SC2*SC2  ‘ ' . 

bC5=SC4*SC 

....  PlD2=PI/2. 


m WMlABli  copy 
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3b* 

*7* 

jd* 

d9* 

Spi02=SQRT  (PU/2) 

VAR=CfM*PID2*SPI02*3,00*S 

SJGSSQRT(VAR) 

DO  26  NTS»2=1»NP1  

40* 

41*  

TS2=NTS2*,1D0+SPMT 
TS2R=TS2*CF 

42* 

TIR=TI*CF 

__43* 

H1R=HI*CF 

44* 

TS1R=TS1*CF 

...  4b*  ; 

hsir=hsi*cf 

46* 

HS2R=HS2*CF 

47* 

AI=S1N(T1R)*C0S(HIR) 

40* 

bl=SlN(UR)*SXN(HIR) 

49* 

Cj=COS(TIR) 

50* 

ASl=SlN(TSlR)*COS(HSiR) 

51* 

bSl=SlN(TSlR)*S!N(HSlR) 

t>2* 

CS1=C0S(TS1R) 

63* 

As2=b IN ( TS2R ) *COS ( HS2R ) 

59* 

6S2=S>  I N (1 S2R ) * S I N ( HS2« ) 

65*  Cs2=C0S(TS2R) 

tjb*  Al=AI+ASi 

57* «JL=bX-*-BSl  - - 

t>d*  Cl=Ci+CSl 

59*  A2=AI+AS2 

60*  B2=6l4-bS2 

61*  C2sCl-fCS2 

62*  aS=(A1+A2)/2# 

_63*  _ . BSS(Wl+B2)/2*  

64*  ArGA1=WN*A1*PL 

6t>*  ARGA2=kN«A2*PL 

6b*  ARGbl=WN*Bl*PL 

67*  AKGB2=WN*B2*PL 

06*  ARGDA=WN*PL*(A1-A2) 

. 69* ARGBbsWN*PL*(Bl-B2) 

70*  AR0A2=ARgDA*AKG0A 
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ArDB2=ARgDB*AKGDB 
SAUtXP  <-ARGAl*ARGAl/2. ) 

SBl=tXP(-ARGBl*ARGBl/2.) . ..  . . 

SA2=tXP ( -ARGA2*ARGA2/2 • ) 

5>B2=EXP  ( -ARGB2*ARGB2/2 • ) 

Cc=SA1*Sa2*SU1*SB2 

CcsCC*GNP  . - . 

ia  omeg=sqrt(g*«<n)*f‘*(Ai»bi) 

IF  (OhEG)  19»20»  19 1 

19  AlOM=CFM*EXPl-2,*G*G/<OMfcG*OMEG*S2) ) 

_ Al0M=Al0M/(0MtG**6)  

CiNF=SQR| ( WN*6) * (PI**3) *PL2*2, 

C IN=C INF*C 1*C2*EXP ( -ARDA*/4 , ) *EXP (-ARDB2/4* ) 
0M1 2=SGR  T(  G*«N*  )*Fl M AS  * BS> 

Om122=OM12**2  . _ 

A12M=CFM*EXP ( -2, *G*G/ < 0Mi22*S2 ) ) 

A12M=A12M/(0M12**6) 

Fl2=F4(AS»BS)**3 

GO  TO  21  ______  

20  ClN=0, 

21  SA12=SA1*SA1  „ 

SB12=SB1*SB1 

_ SA22=SA2*SA2 

Sfc>22=SB2*SB2 

CH=SA12*SB12 

Cll=CH*GNF 

C22=SA22»SB22  _ _ ! 

C22SC22*«NF 

_ F11=(F4<A1»B1))**3  __  

lF(P4(AlrBl))51»50»5i 

51  CI11=CINF*C1*C1*A10M/F11 
GO  TO  52 

50  CI11=0,  _ 

52  0NiEG2=SQKT(WN*G)*F4(A2»B2) 

_ IF  ( ONiEGg ) 23  ’ 24  » 23 4 . 
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liH* 

lib* 

lib* 

H7* 

lib* 

U9* 

1*0* 

1*1* 

1*2* 

1*3* 

1*4* 

1*5* 

l*b* 

1*7* 

l*b* 


lgb*  23  0MEG*2=0MEG2*0MEG2 

lg7*  A20M=CFM*EXP(-2.*G*G/(0MdG22*S2) ) 

- lgb*  A20M=A20M/(0MEG2**6) 

ly9*  F22=(F4(A2»52))**3  _ 

U0»  . CIN=CIN*A1*M/F12 

HI* GXY=CC+CIN  

U2*  C |22=C  INF  *C2*C2*A20M/F22 

U3*  _ _ 60  TO  25  

24  CI22S0, 

25  6XX=Cll+cIll 1 

6YY=C*24-cI22 

C0EN=SQRT(6XX)*S6RT(gYY)  

COHS&XY/CDEN 

C0(NTS2)=ABS(C0H) 

BET=4N*C1*SI6 

__  BfcT12=V»N**Ci*C2*VAR  

*KITE(4'43)SK'SIG’BET'8EU2'TS2'C0H 

26  CONTINUE  . ...  _ 

00  40  NTS2=1’NP1 

Xi  ( NTS2 ) sNTS2* • loO+SpMT  ....  . 

Y1(NTS2)=C0<NTS2) 

40  CONTINUE 

Xl (Nl)sSP 


1*9*  ..  xi(n*)=2.  _____  : - 

150*  Y1(N1)=0. 

lwl*  . Yl(N2)s,2  . ■ — --- 

132*  CALL  LINE(Xl'Yl'NT'l'O'O) 

1^3*  . 0©  TO  30 — . 

lo4*  32  CALC  AXI5(0*»0*  »12HThETA  IN  0E6#-12# 4, #0. *X1 <N1 ) #X1 (N2> » 10, ) . 

13b*  C ALL  AX  I b < 0 . » 0 , » 9HC0hEREnCE  * *9  # 05  • » 90  . , Y1 1 Ni ) , y 1 < N2 ) 1 1 0 , ) 

13b*  CALL  PLOTC20.»0.#-3) 

157*  _ CALL  PLOT  t0.»0,»993) . 

13b*  STOP 

159*  _ ENO -r 

• t l 


